Accepted to the ApJ 

Preprint typeset using LAT^X style cmulatcapj v. 11/10/09 



THE STAR-FORMATION HISTORIES OF Z ~ 2 DOGS AND SMGS 

R. S. Bussmann 1 ' 2 , Arjun Dey 3 , L. Armus 4 , M. J. I. Brown 5 , V. Desai 4 , A. H. Gonzalez , B. T. Jannuzi 3 , J. 

Melbourne 7 , B. T. Soifer 4 ' 7 
Accepted to the ApJ 

ABSTRACT 

The Spitzer Space Telescope has identified a population of ultra-luminous infrared galaxies (ULIRGs) 
at z ~ 2 that may play an important role in the evolution of massive galaxies. We measure the stellar 
masses (M*) of two populations of Spifaer-selected ULIRGs that have extremely red R — [24] colors 
(dust-obscured galaxies, or DOGs) and compare our results with sub-millimeter selected galaxies 
(SMGs). One set of 39 DOGs have a local maximum in their mid- infrared (mid-IR) spectral energy 
distribution (SED) at rest-frame 1.6/zm associated with stellar emission ("bump DOGs"), while the 
other set of 51 DOGs have power-law mid-IR SEDs that are typical of obscured AGN ("power-law 
DOGs"). We measure M» by applying Chariot & Bruzual stellar population synthesis models to 
broad-band photometry in the rest-frame ultra-violet, optical, and near-infrared of each of these 
populations. Assuming a simple stellar population and a Chabrier initial mass function (IMF), we 
find that power-law DOGs and bump DOGs are on average a factor of 2 and 1.5 larger than SMGs, 
respectively (median and inter-quartile M* values for SMGs, bump DOGs and power-law DOGs are 
log(M»/M Q ) = 10.42j;°;3g, 10.62±°;32, and lOJlT^;^, respectively). More realistic star-formation 
histories drawn from two competing theories for the nature of ULIRGs at z ~ 2 (major merger 
vs. smooth accretion) can increase these mass estimates by up to 0.5 dex. A comparison of our 
stellar masses with the instantaneous star-formation rate (SFR) in these z ~ 2 ULIRGs provides a 
preliminary indication supporting high SFRs for a given , a situation that arises more naturally in 
major mergers than in smooth accretion powered systems. 

Subject headings: galaxies: evolution — galaxies: fundamental parameters — galaxies: high-redshift 
— submillimeter 



1. INTRODUCTION 

Ultra-luminous infrared galaxies (ULIRGs) are de- 
fined to have extremely high infrared (IR) luminosities 
(Lm > 10 12 Lq). These luminosities require signifi- 
cant dust heating, usually thought to arise from extreme 
episodes of star-formation (M > 100 Mq yr^ 1 ) or accre- 
tion onto super-massive black holes. These objects are 
rare in the local Universe, yet they have been associated 
with a cr itical phase of galaxy evolution linking merg- 
ers (e.g., lArmus et all 119871; iMurphv et"aH 11996ft with 
quasar s and red, dead elliptical galaxies ([Sanders et al.l 
ll988allbT ). ULIRGs are more commonplace in the dis- 
tant Universe, to the extent that they contribute a sig- 
nificant component of the bolo metric luminosity density 
of the Universe at z > 1 (e.g. iFranceschini et alJ l2001t 
iLe Floc'h et~aLll2005t iPerez-Gonzalez et aLl^OoHT . This 
realization implies that ULIRGs may represent an impor- 
tant evolutionary phase in the assembly history of mas- 
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sive galaxies and has inspired a host of new techniques 
for identifying ULIRGs at z > 1. 

The two most successful techniques for identifying 
high-redshift ULIRGs rely on selection at either mid- 
infrared or far-infrared wavelengths. Surveys at 24/im 
with th e Multiband Ima ging Photometer for Spitzer 
(MIPS; iRieke et al.l 120041) instrument for the Spitzer 
Space Telescope have been remarkab ly successful for 
the mid-IR identification of ULIR Gs (lYan et alJ 120041: 
Houck et all 120051: IWeedman et all I2006bb iFi ore et al.1 
20081: iDev et all 120081 : iFiore et all 12009( 1. In particular, 
Dev et all (j2008f > select sources from the 9 deg 2 NOAO 
Deep Wide-Field Survey (NDWFS) Bootes field that sat- 
isfy R ~ [24] > 14 (Vega magnitudes; ~F 2 4 t _ tm /F R > 
1000) and F 2 4 Mm > 0.3 mJy These objects are called 
dust-obscured galaxies (DOGs), lie at z ~ 2 ± 0.5 
dHouck et alJ 120051: IWeedman etall l2006al iDesai et al.1 
120091 Soi fer et al., in prep. , 2011) have ULIRG luminosi- 
ties (e.g. iBussmann et al.l[2009b( l. have a space de nsity 
of (2.82 ±0.05) x 10" 5 h? m Mpc" 3 l|Dev et al.l 1 2 (M l. and 
inhabit dark matter haloes of mass Mdm ~ 10 12 3 Mq 
(|Brodwin et all 120081) . These results show that DOGs 
are undergoing a very luminous, likely short-lived phase 
of activity associated with the growth of the most mas- 
sive galaxies. 

In addition, DOGs can be divided into two groups ac- 
cording to the nature of their mid-IR spectral energy dis- 
tribution (SED): those with a peak or bump at rest- frame 
1.6/im, likely produced by the photospheres of old stars 
("bump DOGs"), and those dominated by a power-law 
in the mid-IR ("power-law DOGs"). The SED shapes, 
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as w e ll as spectroscopy in the near- IR ([Brand et al 



20071 ISaiina et all 120081) and m id-IR (lYan et al 



Saii na et all 12007): iFarrah et all 12001 iDesai et al.l 



2007 



2009; 



Huang et al.ll2009f ) indicate that the bolomctric luminosi- 



ties of bump DOGs are dominated by star-formation, 
while those of power-law DOGs are dominated by ob- 
scured active galactic nuclei (AGN). This implies that 
the phase of DOG activity is characterized by both vig- 
orous stellar bulge and nuclear black hole growth. 

Another method of selecting high redshift ULIRGs is 
imaging at sub-millimeter (sub-mm) wavelengths. The 
advent of th e Sub-mm Common User Bolometer Ar- 
ray (SCUBA: IBolland et aJJfl999h has allowed wide-field 
surveys at 850/^m which have identified hundreds of 
sub-millimeter selected galaxies (SMGs). These objects 
have similar redshifts (z = 2.2 ± 0.5), number densities 
(n ~ 9 x 10~ 6 h% Mpc" 3 ; iChapman et~aT1 [20051), and 
clust ering properties (Mdm ~ 10 12 2 M©; iBlain et ah! 
l200l as DOGs. 

The fact that SMGs and DOGs have similar prop- 
erties suggests they might be related in an evolution- 
ary sequence simi lar to that of ULIR Gs in the lo- 
cal Universe (e.g. Sanders et al.l 1198 8a). It has been 



hypothesized that such a sequence does indeed exist 
(IDev fc The NDWFS/MIPS Collaboration I2009T) . and 
that DOGs function as an important intermediate stage 
between gas-rich major mergers and quasars at z ~ 
2 (which have similar clus tering properties as DOGs 
and SMGs; iHopkins et al.l [2005 IBrodwin et all l25ol 
iShen et al.ll2009T ) . One intriguing potential piece of sup- 
port for this idea comes from measurements of Ha line 
strengths, which indicate that power-law DOGs have 
lower star-formation rates (SF Rs) by an order of mag- 
nitude compared to SMGs ([Melbourne et al1l2011| ). 

A theoretical understanding of how this evolution- 
ary sequence might occur has recently been advanced 
using iV-body/smoothed particle hydrodynamic simula- 
tions combined with 3D polychromati c dust radiative 
transfer models ([Narayanan et al1l2010l) . In these mod- 
els, simulations are used to follow the evolution of the 
SED of both isolated disk galaxies and major merg- 
ers. These authors find that simulated systems with 
-p24/im > 0.3 mJy are associated with gas-rich (/ g « 0.4) 
major mergers with a minimum total baryonic mass of 
Mb ~ 3 x 10 11 Mq. While there is significant varia- 
tion associated with different viewing angles, initial or- 
bital configurations, etc., the typical simulated major 
merger achieves peak SFRs of ~ 1000 M© yr _1 at the 
beginning of final coalescence when tidal torques funnel 
large quantities of gas in to the nucleus of the system 
(jMihos fc Hernq uist 1996). This period is also when the 
system is brightest at sub-mm wavelengths and thus can 
be selected as an SMG. 

At the same time, central inflows begin to fuel the 
growth of a supermassive black hole. Approximately 
100 Myr after the peak SFR, the black hole accretion 
rate peaks (at about 1-2 M© yr -1 ). The simulations 
include a prescription for active galactic nucleus (AGN) 
feedback that helps terminate star-formation (along with 
consumption of the gas by star- formation). In these mod- 
els, this period of AGN feedback coincides with the DOG 
phase (F 2 4^ m /Fji > 1000). As the gas and dust are con- 
sumed by star-formation, optical sightlines open up and 



the system can be optically visible as a quasar. The evo- 
lutionary progression in the simulations is driven by ma- 
jor mergers and proceeds from SMG to DOG to quasar 
to red, dead, elliptical galaxy (illustrated qualitatively in 
the top panel of Figure [TJ . 

Alternative theories for the formation of SMGs which 
do not i nvolve major mer gers have also been advanced 
recently ([Dave et al.ll2010D . These studies rely on numer- 
ical simulations of cosmological volumes and select SMGs 
as the most actively star-forming systems that match 
the observed number densities of SMGs. The objects in 
the simulations that are designated as SMGs have stel- 
lar masses in the range M* = (1 — 5) x 10 11 M© and 
SFRs in the range 200-500 M© yr -1 . These SFRs are 
a factor of 3 low er than wha t is inf erred observationally 
in SMGs, which iDave et"al"1 ([20101 ) attribute primarily 
to systematic effects in the SFR calibration (in particu- 
lar, a "bottom-light" initial mass function requires lower 
SFRs to produce the observed IR luminosities of SMGs) . 
Because the star-formation histories (SFHs) which pro- 
duce these simulated SMGs do not involve major merg- 
ers, they are referred to here as "smooth accretion" SFHs 
(a qualitative illustration of this SFH is given in the bot- 
tom panel of Figure [1} . 

Studies attempting to connect the mid-IR and far- 
IR selected ULIRG population at high redshift have 
so far focused on th eir basic properties such as bolo- 
metric luminosities (ISaiina et al.l 120081 iCoppin et all 



2008t lLonsdale et all 120091: iBussmann et al.) I2009bl; 



Fiolet et all 12009ft clust ering strengths ([Blain etldl 



20041 : IBrodwin et al.ll2008l ). and mor phologies. In partic 



ular, high-spatial resolu tion imaging (Das vra et alll 2008: 
Melbourne et all 120081 l2009t Bussmann et al.l 12 009a: 
Swinb ank et al.l I2010D and dyna mics ([Tacconi et all 
20061 120081 : Melbourne et al.l 1201 ID have shown no dis- 
tinction in axial ratio that might be suggestive of ori- 
entation effects; instead, these studies have identified 
morphological trends which are consistent with an evo- 
lutionary scenario driven by major mergers in which 
sources that show a bump in their mid-IR SED (i.e., 
bump DOGs and most SMGs) evolve into those with 
a po wer-law dominated mi d-IR SED (i.e., power- law 
DOGs lBussmann et al.l[20TT[ ). To test the origins of these 
sources further, it is imperative to use alternative, com- 
plementary methods of constraining the SFHs of DOGs 
and SMGs at z ~ 2. 

This paper is focused on one such technique: stellar 
population synthesis (SPS) modeling of broad band pho- 
tometry of DOGs and SMGs with known spectroscopic 
redshifts. The primary goal of this study is to place 
the tightest constraints possible given the existing data 
on the stellar masses (M„) and SFHs of bump DOGs, 
power-law DOGs, and SMGs using a uniform SPS mod- 
eling analysis with common model assumptions and fit- 
ting techniques for each ULIRG population. There are 
several reasons to pursue this goal. 

First, constraints on the Af* values and SFHs of 
Spitzer-selected ULIRGs are limited to a few studies 
that have focused on bump sources ([Berta et all [2007; 
lLonsdale et all [2009D . In contrast, the constraints on 
M* and SFHs presented here for power-law DOGs are 
the first such results for this potentially very important 
population of galaxies. If power-law DOGs do not have 
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Figu re 1. Cartoon picture illustrating tw o possible evolutionary paths for massive galaxies at 



2 (adapted from 



Dcy & The NDWFS/MIPS Collaboration 2009). Top: (1) A gas rich major merger leads to an intense, dust-enshrouded phase of star- 
formation. (2) Energetic feedback, possibly from the growth of a central super-massive black hole, heats the dust and gas, cutting off 
star-formation. (3) Depending on the relative timescales of AGN fuelling, dust dissipation, and star formation, the system may be briefly 
visible as a quasar before settling on the red sequence. Bottom: An alternative scenario in which massive galaxies are assembled via smooth 
accretion of gas and small satellites along filamentary structures (some mechanism is still needed to quench star-formation; in this cartoon 
picture, steps (2) and (3) are assumed to be the same as in the major merger driven scenario). One goal of this paper is to test the two 
different possibilities illustrated in step (1) of this diagram using the stellar masses and star-formation rates of high redshift ULIRGs. 



significantly different masses than SMGs or bump DOGs, 
this might imply that the power-law phase occurs during 
the same time that most of the mass in stars is being 
built up. If the power-law is a signature of black hole 
growth, then this would mean that the stellar mass and 
black hole mass are likely being assembled during the 
same period of dust-obscured, intense star-formation. A 
uniform analysis of all three populations is necessary to 
test this hypothesis. 

Second, while SPS modeling methods have become 
more sophisticated, stellar mass results for a given pop- 
ulation haye not necessarily converged. For example, 
iBorvs et al.1 (120051) use Sp itzer/Infr&red Array Camera 
(IRAC: iFazio et al.1 l200l data to infer average SMG 
stellar masses of M* » 2.5 x 10 11 M o re re - 

cently, iDve et all (|2008t l and IMichalowski et all (|2010f l 
have found median stellar masses for SMGs of M* = 6.3 x 
10 n M m and 3.5 x 10 n M , respectively. A new study 
bv lHainline et all (1201 ID using essentially the same data 
set as IMichalowski et al.l (|2010l ) finds significantly lower 
median SMG stellar masses of M* = (7 ± 3) x 10 10 M . 
Finally, measurements of the width of CO emission lines 
in 12 ULIRGs at z ~ 2 have provided a median dynam- 



ical m ass estimate of Md yn ~ 2 x 10 11 M (jEngel et al.l 
l20inh . These sources typically have high gas fractions 
of « 0.5, implying that the stellar masses should be 
M* < 10 11 Mq. This emphasizes the significant sys- 
tematics that affect stellar mass estimates based on SPS 
modeling and underscores the need for a uniform analysis 
when comparing different ULIRG populations. 

Third, the disagreement in observed stellar masses has 
significant bearing on theoretical models for the forma- 
tion of high redshift ULIRGs. As outlined earlier, the 
cosmo logical hydrodynamical simulations of iDave et al.l 
( 2010) predict that SMGs have large stellar masses that 
are ro ughly consistent with the estima tes of IBorvs et al.l 
420051) and IMichalowski et al.l (120100. but a factor o f 
p» 4 larger than th e estim ates of Mainline" "eTaLl (I20T1 . 
The iHainline et al.l (|2011f ) mass estimates are also some- 
what low er than what is expect ated from merger sim- 
ulations ([Narayanan et all 12010). with the caveat that 
such expectations are highly dependent on the stage of 
the merger, viewing angle, etc. A systematic, uniform 
comparison of the relative stellar mass distributions of 
DOGs and SMGs with simulated SFHs from theoretical 
models for the evolution of massive galaxies represents a 
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significant component of this paper. 

In section [2 we present the data used in this analysis, 
including DOG SEDs from rest-frame ultra-violet (UV) 
to near-IR. Section |3] outlines the general methodology 
and describes the SPS libraries, initial mass functions 
(IMFs), and SFHs that are used in the analysis. We 
present our results in section 01 including constraints on 
stellar masses, visual extinctions, and stellar population 
ages. In section [5j we compare our results with simi- 
lar studies of SMGs and other Spitzer-selected ULIRGs 
and explain the implications of the results for models of 
galaxy evolution. Conclusions are presented in section [6l 

Throughout this paper we assume a cosmology in 
which H =70 km s" 1 Mpc" 1 , D, m = 0.3, and D, A = 0.7. 
All magnitudes are in the AB system. 

2. DATA 

The goal of this paper is to study the relative mass 
distributions of samples of high-z ULIRGs, specifically 
DOGs and SMGs, via population synthesis modeling of 
their rest-frame UV through near-IR SEDs. To minimize 
degeneracies in the models, it is important to limit the 
analysis to sources with spectroscopic redshifts. Thus, 
the present sample consists of ULIRGs with spectro- 
scopic redshifts at z > 1.4 and broad-band photome- 
try from the rest-frame UV through near-IR. The sam- 
ple comprises three main sub-groups: two selected with 
Spitzer at 24/im (DOGs), and one selected with the Sub- 
mm Common User Bolometer Array (SCUBA) at 850/im 
(SMGs). 

2.1. DOGs 
2.1.1. Sample Selection 

For the Spitzer-selected ULIRGs, a total of 2603 DOGs 
satisfying R— [24] > 14 (Vega mag) and F 2 4 Mm > 0.3mJy 
were identified in the 8.6 deg 2 NDWFS Bootes fie l d wit h 
deep Spitzer /MIPS 24^m coverage (|Dev et all 120081 ). 
This paper focuses on the subset of 90 of these objects 
that have known spectroscopic redshifts at z > 1.4 either 
from observations with the Keck telescope (~ 60%, Soifer 
et al., in prep., 2011) or with the InfraRed S pectrometer 
(IRS [Houck et al.l 12004ft on board Spitzer (jHouck et al.l 
l2005t iWccdma net al.ll2006bT) . Spectroscopic redshifts for 
our sample of DOGs are given in Table [TJ 

Figure [2] shows the R— [24] color as a function of 24/im 
magnitude for the subsample studied here (the "spec- 
troscopic sample") in comparison to the overall sam- 
ple of DOGs in Bootes. To optimize the spectroscopic 
detection rate, the spectroscopic sample is biased to- 
wards bright 24/mi sources, although the full range of 
R — [24] colors is sampled. The spectroscopic sample 
consists of 39 star-formation dominated "bump" sources 
(those that show a peak at rest-frame 1.6/im) and 51 
active galactic nucleus (AGN) dominated "power-law" 
sources. Bump and power-law DOGs are separated ac- 
cording to the st atistical criteria given in section 3.1.2 of 
iDev et all (|2008| ). Also shown in this diagram are 53 sub- 
milli meter galaxies (SMGs) with spectroscopic redshifts 
from iChapman et~a l. (2005) (see section |2~2|) . The red- 
shift distributions of these groups of galaxies are shown 
in Figure [3] The positions, R — [24] colors, and nature 
of mid-IR SED for each DOG in the sample are given in 
Table [U 
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Figure 2. R — [24] color vs. 24/xm magnitude distribution for 
DOGs in the NDWFS Bootes field. Gray dots and upward arrows 
show the full sample of DOGs, with and without an JS-band detec- 
tion (2<r limits), respectively. Highlighted are the subsamples with 
spectroscopic redshifts and either a mid-IR power-law SED (PL 
DOGs, red circles) or a mid-IR bump SED (Bump DOGs, blue 
squares). Also shown are SMGs (ora nge stars) with spectroscopic 
redshifts from Chapman et al. (2005) and 24/xm photometry from 
IHainline etall 1120091 '). 

2.1.2. Optical Photometry 

The NOAO Deep Wide Field Survey (NDWFS; 
Uannuzi fc Delll999l ) is a ground-based optical and near- 
IR imaging survey of two 9.3 deg 2 fields, one in Bootes 
and one in Cetus. In this paper, we utilize the optical 
imaging of the Bootes field, conducted using the NOAO 
4m telescope on Kitt Peak. The survey reaches 5<r point- 
source depths in %, R, and I of 27.1, 26.1, and 25.4 
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Figure 3. Redshift distribution of DOGs in the Bootes Field 
with spectroscopic redshifts. The redshift distribution of bump 
DOGs (blue hat ched) is rel atively narrow due to selection effects 
(for details see [Dcsai et al. 2009J), while power-law DOGs (red 
hatched) are weighted towards slightly larger redshifts. Also shown 
is the redshift distrib ution of SMGs (orange filled region) from 
Chapman 'eTaTI piM ). 
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(Vega mag), respectively. The NDWFS astrometry is 
tied to the reference frame defined by stars from the 
United States Naval Observatory A- 2 catalog. NDWFS 
data products are publicly available via the NOAO sci- 
ence archive 0- 

Photometry for each DOG was measured in 4" di- 
ameter apertures, centered on the 3.6/im centroid posi- 
tion meas ured from the Spi tzer Deep Wide -field Survey 
(SDW FS; lAshbv et al.ll2009h imaging data ()Ashbv et al.l 
2009) . Foreground and background objects were removed 
using SExtractor segmentation maps, and the sky level 
was determined using an annulus with an inner diam- 
eter of 6" and a width of 5". The background level 
and photometric uncertainty were computed by measur- 
ing the sigma-clipped mean and RMS of fluxes measured 
in roughly fifty 4" diameter apertures within 1' of the 
target. Aperture corrections were derived using bright, 
non-saturated stars for each of the 27 sub-fields that com- 
prise the NDWFS. 

2.1.3. Near-Infrared Photometry 

The NOAO Extremely Wide Field InfraRed iMager 
(NEWFIRM) has conducted a survey at near-IR wave- 
lengths of the full 9.3 deg 2 Bootes field using the NOAO 
4m telescope on Kitt Peak during the spring semesters 
of 2008 and 2009. The nominal 5er limits of the survey 
within a 3" diameter aperture in J, H, and Ks are 22.05, 
21.3, and 19.8 (Vega mag), respectively. All of the survey 
data are publicly available (Gonzalez et al., in prep.). 

Photometry was computed in the same manner as with 
the NDWFS images (see section [2".1.2l) . Aperture correc- 
tions were computed using bright, non-saturated stars for 
each of the 52 sub-fields that comprise the NEWFIRM 
survey of Bootes. Photometry in the optical and near-IR 
is presented in Table [2] 

2.1.4. Mid-Infrared Photometry 

The SDWFS is a four-epoch survey of roughly 8.5 deg 2 
of the Bootes field of the NDWFS. The first epoch of the 
survey took pla ce in 2004 January as part of the IRAC 
Shallow Survey ()Eisenhardt et aLll2004D . Subsequent vis- 
its to the field as part of the SDWFS program reimaged 
the same area three times to the same depth each time. 
The final co-added images have ha depths (aperture- 
corrected from a 4" diameter aperture) of 19.77, 18.83, 
16.50, and 15.85 (Vega mag) at 3.6yLtm, 4.5/mi, 5.8/mi, 
and 8.0/im, respectively. All SDWFS data are publicly 
available. 

Part of the SDWFS Data Release 1.1 includes band- 
matched catal_c^s_created with Source Extractor (SEx- 
tractor, [Bertin^3rnouti [1996). Astrometry in these 
catalogs is tied to 2MASS positions within 0'.'2. We iden- 
tify DOGs in these catalogs using a 3" search radius, and 
use the values in these catalogs for our flux density mea- 
surements of DOGs. SExtractor underestimates the true 
magnitude uncertainties because it assumes a Gaussian 
noise distribution where noise is uncorrelated. In place 
of the SExtractor-derived values, we determine our own 
estimates of the uncertainty on each flux density mea- 
surement using 4" diameter apertures randomly placed 
within V of each object of interest. Photometry in the 
mid-IR is presented in Table [2j 

8 http://archive.noao.edu/nsa 



2.2. SMGs 

2.2.1. Sample Selection 

For the SCUBA-selected SMGs, we use the sam- 
ple of 53 objects with spectroscopic redshifts at z > 
1.4 (we have removed from the sample three sources 
with extremely blue rest-frame ultra-violet colors as 
well as two sources which were subsequently sh own to 
be spurious det e ctions by lHainline et al.l (|201lD ) from 
IChapman et al.l (|2005| ). These are sources with pre- 
cise positional information derived from Very Large Ar- 
ray 1.4 GHz imaging and redshifts obtained with op- 
tical ground-based spectroscopy with the Keck I tele- 
scope. Their clustering properties indicate they inhabit 
very massive da rk matter haloes (Mdm ~ 10 12 M Q ; 
iBlain et al.l 120041 ) , comparable t o the dark matter halo 
masses of DOGs (|Brodwin et alJl2008t ). 

2.2.2. SMG Photometry 

The broad-band photometry of SMGs used in this pa- 
per has been collected from a variety of sources. B- 
and i?-band photometry wer e obtained with several tele- 
scopes and were presented in IChapman et afl (|2005| ). J-, 
J-, and If -band photometry also wer e obtained w i th sev - 
eral telescopes and were presented in lSmail et al.l {2004). 
These photometry values were derived with 4" diame- 
ter apertures and have been aperture-cor rected. Mid-IR 
photo metry of SMGs was obtained from lHainline et al.l 
(2009), who compute aperture-corrected 4" diameter 
aperture photometry using SExtractor. 

3. STELLAR POPULATION SYNTHESIS MODELS 

Stellar population synthesis (SPS) modeling offers a 
means of constraining the mass and star-formation his- 
tory of a galaxy's stellar population. This section con- 
tains a description of the technique adopted here to ap- 
ply the SPS models to the high-z ULIRG photometry 
outlined in section [2j Additionally, details are provided 
regarding three SFHs and three initial mass functions 
(IMFs) that are used in this paper for testing theories 
for the formation of massive galaxies at high redshift. 
Results from this analysis are presented in section 2] A 
detailed analysis of the differences in M* measurements 
obtained with four SPS libraries may be found in Ap- 
pendix [A] 

3.1. General Methodology 

SPS models are parameterized at minimum by their 
luminosity- weighted age and their stellar mass, M*. The 
attenuation of stellar light by dust adds a third pa- 
rameter, Ay. In all models used here, the simplify- 
ing assumption of a uniform dust screen (Ay ranging 
from to 3) is adopted which obscures the intrinsic 
stellar light according to t he re ddening law for star- 
bursts from iCalzetti et al.l (120001) for wavel engths be- 
tween 0.12 — 2.2 /im and that o f lDraind (2003) for longer 
wavelengths. The available data do not allow constraints 
to be placed on more complex models in which younger 
stars have differ ent dust obscuration prescriptions than 
older stars (e.g.. iCharlot fc Falll[200Tl . 

The broad band photometry used here is not sufficient 
to break the degeneracy between age and Ay (except 
under special assumptions). For this reason, the main 
goal here is to measure the relative M* values of three 
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distinct populations of high redshift ULIRGs (power-law 
DOGs, bump DOGs, and SMGs) using a uniform, self- 
consistent analysis. This will allow the stellar masses 
of these objects to be measured in a relative sense and 
therefore minimize many of the uncertainties discussed 
above (however, note that the masses of the power-law 
DOGs in general are upper limits since the AGN con- 
tribution to the 3.6/im and 4.5/im IRAC channels is un- 
known). Furthermore, competing models of galaxy for- 
mation and evolution make different predictions about 
the stellar mass properties of the most luminous galaxies 
at z ~ 2. The distribution of stellar masses of popu- 
lations of power-law DOGs, bump DOGs, and SMGs is 
therefore (in principle) a viable tool with which to test 
these competing models. 

The approach used here is to apply SPS models 
of varying Ay and age values to generate a pro- 
bility density function for the stellar mass of each 
galaxy, </>(M*, age, Ay). </>(M*, age, Ay) is computed di- 
rectly from the best-fit x 2 value for the given num- 
ber of degrees of freedom, -/Vdof- Since we have 
7 data points and 3 model parameters, -/Vdof = 
4. For a few sources (SST24J 142648.9+332927, 
SMMJ030227.73+000653.5, SMMJ123600. 15+621047.2, 
SMMJ123606.85+621021.4, SMMJ131239. 14+424155. 7, 
SMMJ163631.47+405546.9, SMMJ221735. 15+001537.2, 
SMMJ221804.42+002154.4), no models achieved statisti- 
cally acceptable fits. These systems are assumed to have 
a uniform stellar mass probability density function be- 
tween 10 10 — 10 12 Mq. This has the effect of broadening 
the resulting stellar mass constraints for a given galaxy 
population. Each individual galaxy's 4>(M*, age, Ay) is 
normalized such that it contributes equally to the final 
stellar mass probability density function for that popu- 
lation of galaxies (^pldog, ^BumpDOG, and 0smg)- 

The use of SPS models to determine intrinsic proper- 
ties of galaxies assumes that all of the observed flux is 
emitted by stars. In fact, many of the sources in this 
study have a significant contribution in the rest-frame 
near-IR from obscured AGN (this is especially true for 
the power- law DOGs). Some authors add this compo- 
nent (in the form of a variab le slope power-law) t o their 
SPS modelling efforts (e.g.. Mainline et al.ll20lTI ). Al- 
ternatively, it is possible to minimize the AGN contri- 
bution by considering only the first two IRAC channels 
(i.e., up to observed-frame 4.5/im). We adopt the lat- 
ter approach in this study. For bump DOGs and most 
SMGs, this should provide a reasonably reliable measure- 
ment of the stellar light from these objects. For power- 
law DOGs and those SMGs with power-law tails in the 
near-IR, there still exists a significant possibility that the 
observed-frame 4.5/im light is contaminated by AGN, 
though it should be noted that high-spatial resolution 
imaging with ffST/NICMOS indicates that only 10-20% 
of the rest-frame optic al light is emitted by a p oint source 
m power- law DOGs Jl ussmann et all 12009a] ). For this 
reason, the stellar mass estimates of power-law DOGs 
should be regarded as upper limits on the true stellar 
mass. 

The observed-frame B\y photometry have been ex- 
cluded from the fitting process. These data typically 
probe rest-frame 1500 A and as such are highly sensitive 
to the youngest stellar populations and the detailed ge- 



ometry of the dust distribution surrounding them. The 
most robust model fits were obtained when the Bw pho- 
tometry were not used. 

Only solar metallicity models are tested in this study. 
This is a reasonable assumption, since high-redshift (me- 
dian redshift of 2.4) dusty gal axies have been foun d 
to have near-solar metallicities (Swinba nk et al.l 120041) . 
Moreover, our broad-band SED data do not provide the 
ability to constrain metallicity. The adoption of a sin- 
gle metallicity in SPS modeling typi cally introduces un- 
certainties at the level of 10-20% (|Conrov et al.1 120091 : 
M uzzm et al.H2009f ). which are insignificant compared to 
systematic uncertainties related to the IMF, SFH, and 
age of the stellar population. 

3.2. SPS Star-formation Histories 

One of the most critical adjustable parameters in 
SPS modeling is the star- formation history (SFH). 
iMichalowski et all (|201lD suggest that the use of multiple 
component SFHs, in which different stellar populations 
are allowed to have distinct ages and obscuration, can 
lead to factors of 2-4 difference in best-fit stellar mass. 
The focus in this paper is placed on three distinct SFHs 
that broadly encompass a reasonable range of parameter 
space while maintaining a level of simplicity in accor- 
dance with the quality of the available data. 

The first SFH adopted here is the simplest one possi- 
ble: an infinitely short burst of star-formation at time 
t = during which all the stars of the galaxy are formed, 
followed thereafter by passive evolution. This is called a 
simple stellar population (SSP), and is used commonly 
in SPS modeling in the literature. If the objects un- 
der study here have recently had star-formation shut off 
by some process (e.g., AGN feedback), then the SSP 
model provides constraints on how long ago such an event 
ocurred. Models used here have ages spaced logarithmi- 
cally from 10 Myr up to 1 Gyr. 

The second SFH used in this paper is borrowed from a 
representative simulation of a major merger which under- 
goes a very luminous sub-mm phase (SMG) as well as a 
highly dust-obscured phase (DOG) bef ore star-formation 
is sh ut off by AGN feedback effects fNaravana n et all 
2010). This SFH traces the star-formation rate from 
the beginning of the simulation — before the two gas- 
rich (f g ~ 0.8) disks begin to interact — through the 
period of final coalescence when the SFR peaks near 
1000 Mq yr _1 , to the end of the simulation and a red, 
dead, elliptical galaxy. Models used here have ages 
spaced roughly linearly from 10 Myr to 0.8 Gyr. 

The third SFH adopted in this study comes from cos- 
mological hydrodynamical simulations in which SMGs 
are posited to correspond to the most rapidly star- 
forming s ystems that match the observed number density 
of SMGs (pave et al.ll2010T ). In particular, the SFH and 
metallicity history of the highest SFR simulated SMG 
are used. This object has a SFR of w 150 Af© yr _1 for 
most of the simulation but is boosted to « 500 M Q yr _1 
at z = 2 and reaches a mass of M* = 2.8 x 10 11 M© by 
the same redshift. As nearly all of the mass is assembled 
in a quiescent mode, this SFH is nearly opposite to a 
SSP, in which all stars are formed in a single infinitely 
short burst. Models used here have ages spaced roughly 
linearly over the full range of the SFH , from 10 Myr to 
3 Gyr. Figure [4] shows the SFHs from iNaravanan et al.l 
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Figure 4. Star-formation histories used in stellar population syn- 
thesis models. Dotted line represents high-z ULIRGs identified in 
cosmological hydrodynamical simulations formed via smooth gas 
inflow and accret ion of small satellites (Galaxy A from Figure 4 of 
IDave et al.H2010h . Dashed line represents high-z ULIRGs formed 
via major mergers of two gas-rich disks (Narayanan et al. 2010) 
and has been shifted in time so that the peak star-formation rate 
occurs at 2 s; 2. 



(2010) and IDave etHI (|2010h that are used in this anal- 
ysis. 

3.3. Initial Mass Functions 

Another critical adjustable parameter involved in SPS 
modeling is the IMF. Despite its importance, the detailed 
nature of the IMF in galaxies at high redshift is poorly 
constrained. The relevant parameter space is charac- 
terized here by three different f orms: a Salpet er IMF 
(|Salpeterl [T955h a Chab rier IMF (lChabrierll2003l\. and a 
bottom-light IMF (e.g.. Ivan Dokkuml 120081 : lDavell2008h . 
All of these have a lower mass cutoff of 0.1 Mq and an up- 
per mass cutoff of 100 Mq. The Chabrier IMF has fewer 
low mass stars compared to a Salpeter IMF (and hence a 
lower mass- to- light ratio), while a bottom-light IMF has 
even fewer low-mass stars (and a correspondingly lower 
mass-to- light ratio). 

The contribution of low mass stars to the bottom- 
light IMF is governed by the characteristic mass, m c , 
which controls both the cutoff mass at which the lognor- 
mal form dominates as well as the shape of the lognor- 
mal p art of the IMF itself. In particular, Ivan Do kkum 
(2008) use the color and luminosity evolution of clus- 
ter ellipticals to infer m c ~ 2 M q at z > 4 (however, 
see Ivan Dokkum fc Conrovl [20lol which argues instead 
for a steeper-than-Salpeter IMF slope based upon spec- 
tral features that are strong in stars with M* < 0.3 Mq 
found in elliptical galaxies in the local Universe) . In this 
study, a characteristic mass of m c = 0.4 Mq has been 
adopted, as this value matches both the (very rough) 
estimates for SMGs at z ~ 2 as well as theoretical ex- 
pectations based on a model in which the characteristic 
mass is a function of the CMB temperature: m c oc Tq^ b . 
The effect of such a change in the characteristic mass is 
to produce a Salpeter-like slope at M > 1 Mq and a 
turnover at M ~ 1 Mq. This reduces the number of 
low-mass stars relative to the high-mass ones, thereby 
lowering the mass-to-light ratio relative to the Chabrier 
IMF (for intermediate age stars or younger). 



Since observational constraints on the IMF are not 
readily available, each IMF has been tested with each 
SFH (see section |3~2)) . In the case of the simple stellar 
population (SSP), this provides a measure of the uncer- 
tainty resulting from the unknown IMF. However, for the 
purposes of testing the self-consistency of more compli- 
cated SFHs of ULIRGs at high redshift, it is necessary 
to select certain IMFs for each model. The simu l ations 
of major mergers tested here dNaravanan et al.l 12010ft 
adopt a Kroupa IMF for their radiative transfer, so a 
Chabrier IMF (which is very similar to a Kroupa IMF) is 
what is focused on here. Meanwhile, t he IMF is a free pa - 
rameter in the smooth accretion SFH (Dave "eTaI1l2010lh 
A Chabrier IMF is adopted in this paper for this SFH 
(with an accompanying thorough discussion of the impli- 
cations of a more "bottom-light" IMF), since a Salpeter 
IMF overpredicts the sub-mm fluxes of SMGs. 

4. RESULTS 

This section presents measurements of the stellar 
masses (M*) of bump DOGs, power-law DOGs, and 
SMGs. SEDs for each source may be found in Ap- 
pendix [Bj The nominal fiducial model chosen in this 
paper is the CB07 SPS library with a SSP SFH and 
Chabrier IMF (meaning that we have chosen this as the 
standard by which the other models will be compared), 
and is presented in section |4~T1 In later sections, alterna- 
tive SFHs and IMFs are explored. Although differences 
exist between various SPS libraries in the treatment of 
aspects of stellar atmospheres and evolution, these de- 
tails are sub-dominant to the choice of SFH and IMF 
(for an explanation of this, see Appendix |A~| . For this 
reason, our modeling process does not include marginal- 
ization over SPS library. 

4.1. Simple Stellar Population 

The SSP represents a SFH in which all stars form in 
an infinitely short burst of star-formation and evolve pas- 
sively thereafter. While this is an idealized scenario for 
the formation of massive galaxies, it is worth studying 
since SSPs form the building blocks of more complex 
SFHs and can be used more directly to compare the ef- 
fect of different SPS libraries and IMFs (see section |4~41 
for more details on this last point). 

Figure [5] shows the stellar mass probability density 
function resulting from fitting a SSP (computed with the 
CB07 SPS library and a Chabrier IMF) to each power- 
law DOG, bump DOG, and SMG. All three populations 
have a similar range of acceptable M* values. Power-law 
DOGs tend to be the most massive systems, followed 
by bump DOGs and then SMGs. However, their median 
stellar masses are separated by ~ 0.1 to 0.2 dex, while the 
spread in their distributions are rs 0.3 dex. This implies 
that the differences in stellar mass between the popula- 
tions are suggestive rather than conclusive. Perhaps the 
most interesting feature of this result is that the masses 
of all three populations are not significantly different. 
This may imply that that the power-law phase occurs 
during the same time that most of the mass in stars is 
being built up. If the mid-IR power-law is a signature of 
black hole growth, then this implies that the stellar mass 
and black hole mass are being assembled during the same 
period of dust-obscured, intense star-formation. A low 
mass tail is present in each population which is in fact a 
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Figure 5. Stellar mass probability density function of power- 
law DOGs (red), bump DOGs (blue), and SMGs (orange) de- 
rived using using the CB07 library, a Chabrier IMF, and a sim- 
ple stellar population SFH. The me dian M* value (corr ected 
to a Chabrier IMF) fro m studies by Hainlinc ct al] l|201H ) and 
Michalowski ct al. (2010) are given by the long and short dashed 
lines, respectively. The mass estimates presented here indicate 
both types of DOGs have masses similar to S MGs and are closer 
to the Hainlinc ct al. ( 201 1 ) values than those of Michalowski ct al. 
120101) . 

reflection of the fact that the constraints on the stellar 
mass of a small percentage of each group are weak. The 
median and inter-quartile range of stellar masses for this 
SPS model are given in Table [3l 

One feature of the fitting process that is not shown 
in Figure [5] is the well-known significant degeneracy be- 
tween Ay and stellar age - the broad-band photome- 
try of these high-z ULIRGs can be fit either by young 
(10 Myr) and dusty {Ay ~ 1.5 — 2) stellar popula- 
tions or intermediate age (500 Myr) and less dusty 
(Ay ~ 0.0 — 0.5) stellar populations. Given the large 
quantities of dust that are known to exist in these sys- 
tems based on obs e rvations at longer wavelengths (e.g., 
Kovacs et al.l 120061: iCoppin et al.ll200l IBussmann et alJ 



2009al lLonsdale et all 120091: IKovacs et all 12010^ it is 
unlikely that Ay < 1 solutions are acceptable. In- 
deed, mid-IR spectra of Spitzer-selected ULIRGs gener- 
ally show strong silicate abs orption features in dicative 
of highly obscured sources (jSaiina et al.l l2007t ) . Fur- 
thermore, measurements of Ha and H/3 in a handful of 
sources find strong Balmer decrements implying Ay > 1 
(jBrand et al.ll2007f ). 

Assuming Ay = E(B - V)/R v (where Ry — 3.1) and 
the relation betwee n E(B — V) and the hydrogen column 
density (iV H ) from iBohlin et all (| 19781 ). Ay ~ 1 implies 
iVn ~ 2 x 10 21 cm~ 2 . Under the assumption of a spherical 
shell around the source with radius equal to the effective 
radius (R c g), the dust mass can be estimated from Nn 
using: 



1 



(1) 



where / g( j is the gas-to-dust mass ratio (a ssumed to be 
60, th e value found appropriate for SMGs; IKovacs et al.1 
2006) and fi p is the mean molecular weight of the 



gas (assumed to be 1.6 times the mass of a proton). 
Morphological measurements indicat e these object s have 
typical effective radii of 3-8 kpc ( Dasvra et al.l 120081 : 
IBussmann et al.ll2009ak iDonlev et al.ll2010f ). Alltogether 
this implies Md ~ (0.5 — 3) x 10 8 M Q , depending on 
the si ze of R^ff. I n fact , based on 350/im observa- 
tions, IKovacs et al.l (|2010( ) find dust masses of Md ~ 
(5 - 10) x 10 s M Q for Spitzer-selected ULIRGs with a 
mid-IR bump feature. This suggests that Ay > 1 and 
hence age < 200 Myr models should be preferred. Note 
however that for any given galaxy, we do not have in- 
dependent constraints on Ay and hence have applied no 
priors on this quantity in the fitting process. 

4.2. Merger-Driven Star- Formation History 

One of the major goals of this paper is to go be- 
yond instantaneous burst SFHs (SSPs) and test the self- 
consistency of more complicated SFHs. Two in partic- 
ular th at are tested here are a SFH driven by a major 
merger ((Narayanan et al1l2010l ) and a SFH driven mainly 
by smooth accret ion of gas and nearby small satellites 
(jDave et alJl2010j ). The merger-driven SFH is described 
here, while the smooth accretion SFH is described in sec- 
tion O 

Figure [5] (left panel) shows the stellar mass probabil- 
ity density function for power-law DOGs, bump DOGs, 
and SMGs derived us ing a merger-driven SFH (from 
iNaravanan et all 120101 ) with the CB07 SPS library and 
a Chabrier IMF. The median and inter-quartile range of 
M* values are given for this SFH in Tableland are about 
0.1-0.2 dex larger than the same values derived using a 
SSP and a Chabrier IMF (again the trend in masses is 
that power-law DOGs are the most massive and SMGs 
the least massive, with bump DOGs falling in between). 
Multi-component SFHs in general produce higher mass- 
to-light ratios than SSPs because even a modest amount 
of rest-frame UV emission will strongly constrain the age 
of the SSP to be less than a few hundred million years. 
Such a young stellar population will have a low mass-to- 
light ratio. In contrast, a multi-component SFH can have 
a low mass young stellar component (which reproduces 
the rest-frame UV emission) as well as an old stellar com- 
ponent which boosts the mass-to-light ratio. 

This point is made more clearly in the right panel of 
Figure [H which shows the stellar mass probability den- 
sity function for power-law DOGs, bump DOGs, and 
SMGs derived from the merger-driven SFH but focusing 
on the portion of the SFH when the system is expected 
to be in its ULIRG phase (i.e., maximum SFR). By this 
stage (about 0.7 Gyr into the SFH), the presence of a sig- 
nificant amount of low mass stars increases the inferred 
stellar masses by 0.1-0.3 dex (relative to the masses de- 
rived from the SSP SFH). These mass estimates are also 
reported in Table [3] The increase in our estimates of 
M* is actually mitigated somewhat because the SFR is 
so high that the fraction of very massive stars relative to 
all other stars is higher than at other times in the SFH 
and because we have made the simplest possible assump- 
tion for the dust geometry of a uniform dust screen. In 
reality, the youngest stars should experience greater ex- 
tinction than the older stars. This effect is likely to be 
amplified by the merger, in which the peak SFR occurs 
when all the gas and dust have been dumped into the 
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central, most obscured regions. For this reason, we ex- 
pect that our measurements of M* for the merger SFH 
during this period are likely to underestimate the true 
stellar masses. 

It is somewhat interesting that a bimodal distribution 
in SMG stellar masses appears when one focuses on the 
period of peak SFR in the merger simulation. This bi- 
modality is smoothed out in the left panel of Figure [6l 
which shows the superposition of all ages during the SFH. 
The origin of the bimodality is not entirely clear, but 
is likely due to the presence of a significant number of 
SMGs that are rest-frame UV-bright and therefore are 
found to have relatively low stellar masses. In contrast, 
DOGs are selected to be rest-frame UV-faint and do not 
show this bimodality in stellar masses. 

4.3. Smooth Accretion Star- Formation History 

In the c o smolo gical hydrodynamical simulations of 
IDave et al.l ()2010f ). SMGs are posited to be the max- 
imally star-forming galaxies whose number densities 
match the observed number density of SMGs. This re- 
sults in the typical simulated SMG having a SFH de- 
scribed by a relatively constant SFR of 100-200 M yr _1 
over a period of 3 Gyr and leads to a stellar m ass in these 
system s in the range M© w (1—5) x 10 11 M . IDave et"al"1 
(2010) note that their simulated SFRs are a factor of ~ 3 
lower than the typical values observationally inferred for 
SMGs, and hypoth esize that a "bottom -ligh t" IMF such 
as that proposed bv lvan Dokkum (2001) and|Davl(2008) 
could explain this discrepancy. This type of IMF would 
also have the consequence of modifying the M*/Ly of the 
galaxy, meaning that at a given Ly, the inferred stellar 
mass will be lower than for other IMFs such as Chabrier 
or Salpeter. It is for this reason that the constraints on 
the stellar masses of the high-z ULIRGs with this SFH 
are of particular interest. 

Figure [7] (left panel) shows the stellar mass probabil- 
ity density function for power-law DOGs, bump DOGs, 
and SMGs derived using a SFH driven mainly by smooth 
accretion of gas and nearby satellites (with the CB07 
SPS library and a Chabrier IMF). The median and inter- 
quartile range of -M* estimates are provided in Table [H 
In this case, the median stellar masses of the three pop- 
ulations are separated by « 0.15 dex, with power-law 
DOGs being the most massive and SMGs being the least 
massive (note that this is still well below the typical inter- 
quartile range in the stellar mass estimates of « 0.3 dex). 
In comparison to the SSP SFH, the smooth accretion 
mass estimates are ~ 0.2 dex larger, for similar reasons 
as those outlined at the end of section 14.21 

Restricting the age range of the SFH for the smooth 
accretion model to coincide with the period during which 
the simulated systems are expected to be ULIRGs (i.e., 
at z ~ 2 — 3) leads to inferred stellar masses that are 
larger by 0.3-0.4 dex compared to the SSP SFH (Tabled 
and Figure [3 right panel). As described in section |4~2I 
this is a result of a greater contribution from older stars 
that have higher mass-to-light ratios than younger stars. 

4.4. Variation with IMF 

In SPS modeling, the IMF affects primarily the 
mass-to-light ratio o f the synthesized stellar population. 
iBruzual k Charlotl (12003ft showed that the B — V and 



V — K colors of SPS models distinguished only by their 
IMFs (Chabrier vs. Salpeter) are very similar. On 
the other hand, the Salpeter IMF gives mass-to-light 
ratios that are « 0.2 dex larger than the Chabrier 
IMF. Bottom-ligh t IMFs (such as that advocated by 
van Dokkum 2008) have more complicated mass-to-light 
ratios that depend on both the characteris tic mass (m c ) 
and t he age of the stellar population. Ivan Dokkum! 
(2008) find that for m c = 0.4 M (as adopted here) 
and ages < 1 Gyr, the mass-to-light ratio is lower by 
0.2-0.3 dex compared to a Chabrier IMF. The results of 
this study are consistent with this finding: assuming a 
SSP SFH and this bottom-light IMF, the stellar masses of 
bump DOGs are in the range M* = (0.1-0.6) x 10 11 M©, 
or about 0.3-0.4 dex lower than those inferred using a 
Chabrier IMF. A similar reduction in M* occurs when 
using the bottom-light IMF in conjunction with more 
complicated SFHs such as the merger-driven SFH and 
the smooth accretion SFH detailed in sections 14.21 and 
1431 

5. DISCUSSION 

The focus of this section is to build upon the con- 
straints on the stellar masses and star-formation histo- 
ries of bump DOGs, power-law DOGs, and SMGs pre- 
sented in section [4] Estimates of M* presented here are 
compared with estimates of other dust-obscured high- 
redshift ULIRGs. In addition, implications for models 
of galaxy evolution are presented based upon a compari- 
son of the two theoretical SFHs considered in this study 
(major merger and smooth accretion). 

5.1. Comparing Stellar Mass Estimates of ULIRGs at 

Studies of other Spitzer-selected ULIRGs with a bump 
in the observed-frame mid-IR SED have found median 
stellar masses of M, » 10 11 M w ( f or a Chabrier IMF; 
Berta et al.l 120071: lLonsdale et all 120091: iHuang et all 
20091 ). This is a little more than la higher than the 
median stellar mass found for bump DOGs here. The 
small difference in stellar mass estimates can be fully ac- 
counted for by the choice of star-formation history as 
well as the use in this study of the new CB07 SPS li- 
braries, which have redder near-IR colors and hence tend 
tow ards lower inferred s tellar masses (see section l4~4l and 
also iMuzzin et all 120091 ) . 

Two recent studies of SMGs using stellar population 
synthesis modeling have come to di ffering conclusions 
regard ing their median M*. While IMichalowsk i et al.l 
(|2010ft find a medi an stellar mass of M* » 2 x 10 11 M Q 
(using SEDs from llglesias-Paramo et al.l 120071 and af - 
ter converting to a Chabrier IMF) . lHainline 
find M* = 7 x 10 10 M (ass u ming a Chabrie r IMF and 
models from iMarastonl 12005ft ■ lHainline et~ail (|2011[ ) ar- 
gue that models which do not consider the contribu- 
tion of an obscured AGN in the mid-IR (particularly 
in the 5.8/im and 8.0/im channels of IRAC) can bias 
stellar mass estimates of SMGs upwards by a factor of 
w 2. Our analysis (which excludes these two IRAC 
channels to minimize the contribution from an obscured 
AGN) indicates stellar masses that are closer to those of 
lHainline etlll (|201lD . with median M* = 4.4 x 10 10 M Q 
(for a smooth accretion SFH without constraints on the 
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Figure 6. Similar to Figure [5] but assuming a major merger SFH (see Figure [4] of this paper and Narayanan ct al. 20 101*) . Left: Stellar 
masse s obtained when marginalizing over the full age range of the model. Mass estimates from Hainlinc ct al. (2011) and Michalowski ct al. 
(2010) (dashed and dotted lines, respectively) have been corrected to a Chabrier IMF. Right: Stellar masses obtained when marginalizing 
over only the 2 = 2 period of the model (i.e., the timestep during which the SFR peaks). When only the z = 2 period is considered, the 
inferred stellar masses increase by 0.1-0.3 dex (see Table[3} thanks to the increased contribution from old stars with high mass-to-light 
ratios. The median stellar masses of these z ~ 2 ULIRGs are still about 0.2-0.3 dex lower than expected from the [Narayanan et al. 12010) 
models (solid line). 
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Figure 7. Similar to Figure \E\ but assuming a smooth accretion SFH (see Figure [4] of this paper and IDave eTaLl 120101 ). Left: Stellar 
masses obtained when marginalizing over the full age range of the model. Right: Stellar masses obtained when marginalizing over only the 
z ~ 2 — 3 period of the model. Here, the stellar mass values increase by 0.3-0.4 dex when only the z ~ 2 epoch of the SFH is considered 
due to the increased contribution of low-mass stars whic h have high mass-t o-light ratios. The median stellar masses of these z ~ 2 ULIRGs 
are still about 0.3-0.4 dex lower than expected from the Dave ct al. (2010) models (solid line). 



age of the stellar population, which most closely re- 
sembles the SSP and constant s tar formation histories 
adopted by Mainline et al.l (|2011[ )). Inclusion of the ad- 
ditional two IRAC channels increases our stellar mass 
estimates by 50% (median M* = 6 X 10 10 M ). A simi- 
lar increase is seen when including the 5.8/im and 8.0/im 
channels of IRAC to models which focus on the age of 
the smooth accretion SFH corresponding to the ULIRG 
phase (i.e., z ~ 2). This is somewha t of a smaller effect 
than found by Mainline et al.l (|2011l ). possibly suggest- 
ing that there may be some contamination from AGN 



in the 3.6/Ltm and 4.5/im IRAC channels that we are not 
accounting for as well. 

5.2. Implications for Galaxy Evolution at z ~ 2 

Observational evidence indicates that ULIRGs in 
the local Universe are the product of major mergers 
(jArmus et al.l fl987) and that t hey are connected in an 
evolutionary sense with quasars ()Sanders et al.lfl9"88airbT ) . 
It is tempting to postulate a similar major-merger ori- 
gin for high-redshift ULIRGs. However, conclusive evi- 
dence linking variously selected ULIRG populations to 
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each other and to quasars at high redshift requires 
measurements that challenge current observational ca- 
pabilities. Nevertheless, some tantalizing hints exist 
that suggest these diverse populations are indeed linked. 
First, the clustering strength of DOGs is comparable 
to tha t of both the SMGs and QSOs at similar red- 
shifts (|Brodwin et al.l 12008). Second, the quantitative 
morphologies of DOGs and SMGs are consistent with 
an evolutionary picture in which the SMG phase pre- 
cedes the bump DO G phase, which in tu rn precedes 
the PL DOG phase (jBussmann et al.l 120111 ). However, 
such morphological studies are challenging because of 
surface brightness dimming and dust-obscuration effects, 
which prevent a s traightforward merger identification 
based on imaging (iDasvra et al". 2008; Melbourne et all 
20081 IBussmann et all l2009al Melbourne et al.l 120091 : 
Zamoiski et al.ll201.il IBussmann et al.ll2011| ). 

This study offers an independent means of testing both 
the evolutionary hypothesis as well as the merger hy- 
pothesis via SPS modeling of broad-band imaging in the 
rest-frame UV through near-IR. The approach followed 
in this paper is to test the self-consistency of two dis- 
tinct SFHs. One is characterized by a gas-rich major 
merger which reaches a peak SFR of s» 1000 M@ yr -1 
(|Naravanan et all 12010). The other is characterized by 
smooth accretion of gas and small satellites that typi- 
cally reaches SFR s of « 200 - 300 M Q yr" 1 at z ~ 2 
(pave et al.l[20Toh . 

First, it is worth noting that in the model of 
iNarayanan et a l. (2010), an evolutionary progression ex- 
ists in which SMGs evolve into bump DOGs which evolve 
into power-law DOGs. This process occurs on a short 
time-scale (~ 50 — 100 years), but if it is true then we 
should expect the SMGs to have the youngest stellar 
population, followed by bump DOGs and then power- 
law DOGs. In this case, the relative differences in the 
inferred stellar masses for the three populations become 
more significant and in the expected direction for the 
evolution ary scenario outlin ed above. In comparison, the 
model of iDave et al.l (|2010( ) does not yet include radia- 
tive transfer calculations and so cannot make a prediction 
for an evolutionary scenario between these three popula- 
tions. Because of the short timescales involved in the 
merger simulations and the nature of our seven filter 
broadband photometry dataset, we do not pursue this 
point in a more quantitative manner, but we neverthe- 
less believe that it deserves mentioning. 

Second, the stellar masses reported here are factors of 
2-2.5 lower than expected from the both the merger and 
smooth accretion models tested in this paper. This re- 
flects the large uncertainties inherent in absolute mea- 
surements of stellar mass and indicates that stellar 
masses alone are unlikely to provide a definitive reason 
to fav or either model over the other. Michalowski et al. 
(|2011|) show that the use of multi-component SFHs (i.e., 
multiple stellar populations with varying ages, extinc- 
tions, and masses) in SPS modeling can lead to higher 
inferred total stellar masses by virtue of using the young 
stellar component to match the rest-frame UV flux and 
the old stellar component to match the rest-frame near- 
IR flux. We do not believe the data we have in hand 
(broad-band photometry in seven filters) is sufficient to 
warrant such complex models, but it is nevertheless im- 
portant to recognize that such models are indeed capable 



of implying larger stellar masses than the models we have 
adopted in this paper. 

In addition, the unknown form of the IMF can po- 
tentially insert another factor of 2-4 uncertainty in the 
absolute stellar mass measurements. However, it must 
be emphasized that modifications in the assumed IMF 
will affect not only the inferred stellar masses, but also 
the inferred instantaneous SFRs. Thus, the effect of the 
IMF can be minimized by comparing the stellar masses 
of z ~ 2 ULIRGs to their star-formation rates. Although 
SFRs are not yet well known in Spitzer-selected ULIRGs, 
early evidence indicates that bump so urces may have 
simil a r SFRs as SMGs (~ 1000 M yr" 1 lLonsdale et all 
120091: iKovacs et al.1 12010D . whereas power-law sources 
may have much lowe r SFRs (e.g. ~ 100 M yr , 
Melbourne et al.ll20lil ). 

A galaxy with a mass of M* = 1 x 10 11 Mq and a 
SFR of 1000 M Q yr" 1 has a specific SFR of sSFR = 
1 x 10~ 8 yr -1 . A galaxy with a SFR lower by a factor 
10 will have a sSFR that is also lower by a factor of 10. 
Thus the range for DOGs and SMGs in sSFR is likely 
to be of order 1 — 10 Gyr~ . In comparison, simulations 
of major mergers that produce DOG and SMG behavior 
tend to have sSFR = 6 Gyr -1 . On the other hand, in 
smooth accretion driven simulations, SMGs have sSFR = 
0.7 Gyr -1 . Even if we adopted assumptions regarding 
the SFHs and dust geometry that led to stellar masses 
that were a factor of 2-4 l arger and were thu s consi stent 
with those found by e.g. [Michalows ki et all (|2010f) . the 
range in sSFR values for DOGs and SMGs would still be 
higher than the expectation from the smooth accretion 
model. This is merely a consequence of the fact that 
mergers provide a more ready mechanism to obtain high 
sSFR values than smooth accretion models. 

6. CONCLUSIONS 

In this paper, we analyze the broad-band SEDs of a 
large sample of mid-IR selected (bump and power-law 
DOGs) and far-IR selected (SMGs) ULIRG populations 
with known spectroscopic redshifts and use stellar popu- 
lation synthesis models to estimate self-consistently the 
stellar masses of these three populations. We compare 
our mass estimates with predictions from two competing 
theories for the formation of these systems and examine 
the implications for galaxy evolution. We list our find- 
ings below. 

• The median and inter-quartile range of stellar 
masses for SMGs, bump DOGs and power-law 
DOGs are log(M*/M ) = 10.71±g;|2, 10.62±£g|, 

and 10.42to;3g, respectively, assuming a simple 
stellar population SFH, a Chabrier IMF, and the 
CB07 stellar libraries. The overlap in M* values be- 
tween all three populations is consistent with the 
picture in which they represent a brief but impor- 
tant phase in massive galaxy evolution, with tenta- 
tive evidence supporting a scenario in which SMGs 
evolve into bump DOGs which evolve into power- 
law DOGs. 

• The use of more realistic SFHs in the SPS model- 
ing in which both old and young stars contribute to 
the observed broad-band photometry can increase 
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mass estimates significantly. We show that us- 
ing a major merger driven SFH during its peak 
SFR period (when it is expected to be identi- 
fied as a ULIRG at z ~ 2) leads to median and 
inter-quartile stellar mass estimates for power-law 
DOGs, bump DOGs, and SMGs of log(M*/M ) = 
H-06ioJt 10.88tg;5J, and 10.86±g;^, respectively. 
Using a smooth accretion driven SFH (focusing 
on the predictions at z ~ 2) these values be- 
come log(M*/M ) = 11.20tg;|5, 11.03^;^, and 
11.02±g;|f, respectively. 

• The stellar masses we measure are inconsistent 
with those predicted by both numerical simula- 
tions we have tested (being lower by a factor of 
2-2.5). This indicates that either the simulations 
over-predict the stellar masses of high-z ULIRGs, 
or that one (or more) of the assumptions in our 
SPS models is incorrect. In either case, the stellar 
mass data presented here are by themselves insuf- 
ficient to favor one model over another. However, 
we note that the use of a bottom-light rather than 
a Chabrier IMF may be needed for the SFRs of 
the smooth accretion model to match those that 
are observed. Such a change would decrease our 
mass estimates by a factor of roughly 2 (depend- 
ing on the exact shape of the bottom- light IMF). 
This line of reasoning suggests that, at least for the 
most luminous sources, the smooth accretion model 
has difficulty reproducing the observed far-IR emis- 
sion (i.e., instantaneous SFR) without overestimat- 
ing the observed optical and near-IR emission (i.e., 
stellar mass). 

Estimates of the stellar masses of dust-obscured galax- 
ies at high-redshift are highly dependent on the age of 
the stellar populations within those galaxies. The use of 
multiple component SFHs with different ages can lead 
to significant variations in the inferred stellar mass (e.g., 
iMichalowski et afll2011[ ). We do not consider the broad- 
band photometry in seven filters used here to be suf- 
ficient to explore such complex SFHs. However, in the 
near future, wide-field medium-band photometry surveys 
in the near- IR (e.g., the NEWFIRM Medium-Band Sur- 
vey, NMBS; lvan Dokkum et al.ll2009f) will provide a finer 
sampling of the rest-frame Balmer and 4000 A break and 
significantly improve constraints on the stellar popula- 
tion age in DOGs and SMGs. Further in the future, the 
advent of the James Webb Space Telescope will provide 
high-spatial resolution imaging in the mid-IR and pro- 
vide improved constraints on the amount of stellar emis- 
sion vs. AGN emission in ULIRGs at high redshift. This 
is critical information especially for power-law DOGs, 
but holds significance for bump DOGs and SMGs as well. 

This work is based in part on observations made with 
the Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technol- 
ogy under NASA contract 1407. Spitzer/MIPS guaran- 
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at 24/mi and is critical for the selection of DOGs. 
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APPENDIX 
SPS LIBRARIES 

Four SPS libraries have been tested in t his an alysis of the SEDs of DOGs and SMGs. The first SPS library used in 
this paper is from thelBruzual fc Chariot! ()2003P ) population synthesis li brary. It uses the isq chrone synthesis technique 
(|Charlot fc Bruzual 11991 ) and the Padova 1994 evolutionary tracks (|Girardi e t al. 1996) to compute the spectral 
evolution of stellar populations at ages between 10 5 and 2 x 10 10 yr. The STEllar LIBrary (STELIB iLe Borgne et al.l 
2003) of stellar spectra offer a median resol ving power of 2000 ov er the wavelength range 3200 to 9500 A. Outside this 
wavelength range, the BaSeL 3.1 libraries (|Westera et al.ll2002D are used and offer a median resolving power of 300 
from 91 A to 160/mi. 

The second SPS library used here is an updated version of the lBruzual fc Charlotl (j2003l ) population synthesis library 
(Chariot & Bruzual, private communication, hereafter CB07). The primary improvement included in these models is a 
new prescription for the therm ally pulsing asymptotic g iant branch (TP-AGB) evolution of low- and intermediate-mass 
stars iMarigo fc Girardil (|2007|) and iMarigo et al.l 1)20081 ) . This has the effect of producing significantly redder near-IR 
colors for young and intermediate-age stellar populations, which leads to younger inferred ages and lower inferred 
masses for a given observed near-IR color. These new models otherwise still rely on the Padova 1994 evolutionary 
tracks and the combination of BaSeL 3.1 and STELIB spectral libraries. 

The third SPS l ibrary employed in this pap er is called a Flexible Stellar Population Synthesis library (FSPS; 
IConrov et afll200l l20ia IConrov fc Gundl2^1fll. This library uses the isochrone synthe sis technique as well, but with 
updated evoluti onary tracks (Padov a 2008 IMarigo fc Girardil l2007HMarigo et al.ll2008f ). FSPS adopts the BaSeL 3.1 
spectral library (|Westera et al.|[2002f ) but includes TP- AGB spectra from a compilation of more than 100 optical/near- 
IR spectra spanning the wavelength range 0.5 2.5/mi (|Lancon fc Woodll2000t iLancon fc Mouh cine 2002). One feature 
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of this library that is not available in the others is the ability to i nput a custom IMF (e.g., a "bottom-light" IMF). 

The fourth and final SPS library used here is from [Maraston] (|2005l) . This library adopts the "fuel-consumption" 
approach, in which the integration variable is the amount of hydrogen or helium consumed by nuclear burning during 
a given post-main-sequence phase (unlike the isochrone synthesis approach, in which the integration variable is the 
stellar mass). This library features a strong contribution from TP- AGB stars (» 40% of the b olometric light) for age 
ranges of 0.2 - 2 Gyr. A comparison between this libra ry and that of Bruzual fc Chari ot (2003) found that the near-IR 
colors of z ~ 2 galaxies were better fit by the former ([Maraston et al.ll2006() . highlighting the importance of a proper 
treatment of the TP-AGB phase for intermediate age stellar populations. 

SEDS 

Since every source in this study has a known spectroscopic redshift, it is possible to construct SEDs 
for each source showing the luminosity per unit frequency (L u ) as a function of rest-frame wavelength 
(A r cst)- Figures [H HI UHl show the SEDs for power-law DOGs, bump DOGs, and SMGs, respectively. 
Also shown in this diagram is the best-fit synthesized stellar population model (CB07, simple stellar pop- 
ulation, Chabrier IMF). Inset in each diagram is the stellar mass probability density function. In a few 
cases (SST24J 142648.9+332927, SMMJ030227.73+000653.5, SMMJ123600. 15+621047.2, SMMJ123606.85+621021.4, 
SMMJ131239. 14+424155. 7, SMMJ163631.47+405546.9, SMMJ221735.15+001537.2, SMMJ221804.42+002154.4), no 
acceptable model was found within the probed region of parameter space. In the subsequent analysis, these systems 
are assumed to have a uniform stellar mass probability density function between 10 10 — 10 12 M©. 

Power-law DOGs have the brightest rest-frame near-IR luminosities, with luminosities at 3/im approaching vL v = 
10 12 Lq. This represents a near-IR excess of a factor of 3-5 compared to bump DOGs a nd SMGs. S uch an excess 
is an indicator of thermal emission from an obscured nuclear source (i.e., obscured AGN; iRiekd 11978ft . Meanwhile, 
bump DOGs and SMGs have rest-frame optical and near-IR SEDs that qualitatively match the shape of the synthezied 
stellar population shown in Figures [9] and [TQl This is consistent with the notion that this part of the SED of these 
objects is dominated by stellar light. 

Relative to their rest-frame near-IR luminosities, SMGs show a rest-frame UV excess compared to bump DOGs and 
power-law DOGs. This is likely the result of a selection effect, but the physical implications are unclear. Possible 
explanations include a difference in dust obscuration or in the luminosity weighted-age of the stellar population. 
Resolving this issue may require deep, high spatial resolution imaging of SMGs in the rest-frame UV, optical, and 
near-IR (currently, only UV and optical imaging is available and only for a handful of sources; e.g. iConselice et al.l 
[2001 ISwinbank et al.ll20Toh . 
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Table 1 

Basic DOG Spectroscopic Sample Properties 



ID R.A. (,12000) Dec. (J2000) z a Bump/Powcr-law R - [24] 



SST24 


,1142538 


.2 + 351855 


216 


.4089050 


35 


.3156586 


2.26 


Pt 




r-law 


> 15.6 


SST24 


J142541 


.3 + 342420 


216 


.4219513 


34 


.4056931 


2.194 


p, 




r-law 


14.7 


SST24 


J142554 


.9+341820 


216 


.4792328 


34 


.3057480 


4.412 


Pt 




r-law 




SST24 


,1142607 


.8 + 330425 


216 


.5326385 


33 


.0739212 


2.092 


p, 




r-law 


14. 1 


SST24 


,1142622 


.0+345249 


216 


.5918884 


34 


.8804398 


2.00 




Bu. 




15.0 


SST24 


.1142626 


.4+344731 


216 


.6102295 


34 


.7919617 


2.13 


Pt 




r-law 


> 15.7 


SST24 


,1142637 


.3 + 333025 


216 


.6558075 


33 


.5071220 


3.200 


p, 




r-law 


> 14.7 


SST24 


,1142644 


.3 + 333051 


216 


.6846313 


33 


.5143967 


3.312 


Pt 




r-law 


14.9 


SST24 


,1142645 


.7+351901 


216 


.6904144 


35 


.3169899 


1.75 


p, 




r-law 


> 16.3 


SST24 


J142648 


.9+332927 


216 


.7039337 


33 


.4908333 


2.00 


Pt 




r-law 


15.7 


SST24 


,1142652 


.5 + 345506 


216 


.7188568 


34 


.9181824 


1.91 




Bu: 




15.0 


SST24 


J142653 


.2 + 330221 


216 


.7218781 


33 


.0391388 


1.86 


Pt 




r-law 


15.8 


SST24 


J142724 


.9+350824 


216 


.8541260 


35 


.1399765 


1.70 




Bu: 




> 14.8 


SST24 


J142748 


.4+344851 


216 


.9518738 


34 


.8142471 


2.200 


Pt 




r-law 


14.6 


SST24 


J142759 


.8+351243 


216 


.9991150 


35 


.2118530 


2.100 


p, 




r-law 


> 15.4 


SST24 


,1142800 


.6 + 350455 


217 


.0028992 


35 


.0819473 


2.223 


P. 




r-law 


14.7 


SST24 


J142804 


.1 + 332135 


217 


.0172119 


33 


.3596916 


2.34 




Bu. 




> 15.8 


SST24 


J142810 


.5 + 352509 


217 


.0439453 


35 


.4192238 


1.845 


P. 




r-law 


14.8 


SST24 


J142814 


.2 + 352245 


217 


.0593109 


35 


.3795052 


2.387 


Pt 




r-law 


14.2 


SST24 


J142815 


.4+324720 


217 


.0640869 


32 


.7887993 


2.021 


P, 




r-law 


15.1 


SST24 


,1142827 


.9+334550 


217 


.1163635 


33 


.7639198 


2.772 


Pt 




r-law 


15.4 


SST24 


,1142832 


.4+340849 


217 


.1351166 


34 


.1473694 


1.84 




Bu. 




13.8 


SST24 


J142842 


.9+342409 


217 


.1790771 


34 


.4030418 


2.180 


P, 




r-law 


15.1 


SST24 


,1142846 


.6+352701 


217 


.1942139 


35 


.4504471 


1.727 




Bu: 




> 15.3 


SST24 


J142901 


.5 + 353016 


217 


.2565460 


35 


.5044174 


1.789 


P, 




r-law 


> 14.7 


SST24 


J142920 


.1 + 333023 


217 


.3341827 


33 


.5063858 


2.02 




Bu: 




14.0 


SST24 


J142924 


.8 + 353320 


217 


.3533783 


35 


.5559425 


2.73 


P, 




r-law 


> 15.9 


SST24 


,1142928 


.5 + 350841 


217 


.3685455 


35 


.1448898 


1.855 




Bu: 




> 14.4 


SST24 


J142931 


.3+321828 


217 


.3808136 


32 


.3076057 


2.20 


P. 




r-law 


> 15.7 


SST24 


J142934 


.2 + 322213 


217 


.3932343 


32 


.3701096 


2.278 


Pt 




r-law 


15.2 


SST24 


J142941 


.0+340915 


217 


.4209595 


34 


.1542397 


1.90 




Bu: 




> 14.7 


SST24 


,1142951 


.1 + 342042 


217 


.4629822 


34 


.3447685 


1.77 




Bu 


'lip 


> 14.7 


SST24 


,1142958 


.3 + 322615 


217 


.4930878 


32 


.4376068 


2.64 


p, 




r-law 


15.6 


SST24 


J143001 


.9+334538 


217 


.5076904 


33 


.7603149 


2.46 


Pt 




r-law 


16.2 


SST24 


,1143020 


.4+330344 


217 


.5855865 


33 


.0622444 


1.482 




Bu: 




> 15.2 


SST24 


J143022 


.5 + 330029 


217 


.5941925 


33 


.0080185 


3.15 


Pt 




r-faw 


> 15.7 


SST24 


J143025 


.7+342957 


217 


.6072998 


34 


.4992828 


2.545 


p, 




r-law 


15.4 


SST24 


,1143028 


.5 + 343221 


217 


.6188049 


34 


.5392456 


2.178 


P. 




r-law 


15.1 


SST24 


J143102 


.2 + 325152 


217 


.7593689 


32 


.8645210 


2.00 


P> 




r-law 


> 15.8 


SST24 


J143109 


.7+342802 


217 


.7908020 


34 


.4673615 


2.10 


P. 




r-law 


15.7 


SST24 


J143135 


.2 + 325456 


217 


.8971863 


32 


.9158325 


1.48 


Pt 




r-law 


14.7 


SST24 


J143137 


.1 + 334501 


217 


.9042053 


33 


.7503319 


1.77 




Bu: 




14.8 


SST24 


J143152 


.3 + 350030 


217 


.9683838 


35 


.0082169 


1.52 




Bu. 


'lip 


14.6 


SST24 


J143201 


.8 + 340408 


218 


.0076141 


34 


.0688477 


1.857 


p, 




r-law 


14.5 


SST24 


J143216 


.8 + 335231 


218 


.0702515 


33 


.8754730 


1.76 




Bu 




> 14.8 


SST24 


,1143225 


.3 + 334716 


218 


.1057739 


33 


.7878914 


2.00 


p, 




r-law 


> 15.9 


SST24 


J143242 


.5 + 342232 


218 


.1771698 


34 


.3757019 


2.16 


Pt 




r-law 


> 15.5 


SST24 


J143251 


.8 + 333536 


218 


.2159729 


33 


.5932732 


1.78 


p, 




r-law 


> 15.3 


SST24 


J143312 


.7+342011 


218 


.3028564 


34 


.3364716 


2.119 


P. 




r-law 


15.3 


SST24 


J143315 


.1 + 335628 


218 


.3133240 


33 


.9411583 


1.766 


P> 




r-law 


14.2 


SST24 


J143318 


.8+332203 


218 


.3284149 


33 


.3674889 


2.175 


P. 




r-law 


14.6 


SST24 


J143321 


.8 + 342502 


218 


.3410492 


34 


.4173508 


2.09 




Bu. 




14.2 


SST24 


J143324 


.3 + 334239 


218 


.3508911 


33 


.7109337 


1.93 




Bu. 


mp 


14.3 


SST24 


J143325 


.8 + 333736 


218 


.3575897 


33 


.6268959 


1.90 


Pt 




r-law 


15.4 


SST24 


,1143330 


.0+342234 


218 


.3752289 


34 


.3762436 


2.082 


p, 




r-law 


15.2 


SST24 


J143331 


.9+352027 


218 


.3831787 


35 


.3409195 


1.92 




Bu 




14.3 


SST24 


J143332 


.5 + 332230 


218 


.3855133 


33 


.3750801 


2.778 




Bu. 


mp 


> 15.3 


SST24 


,1143335 


.9+334716 


218 


.3996735 


33 


.7877769 


2.355 


Pt 




r-law 


14.5 


SST24 


,1143349 


.5 + 334601 


218 


.4567871 


33 


.7671394 


1.87 




Bu 




> 14.7 


SST24 


J143353 


.7+343155 


218 


.4738007 


34 


.5321503 


1.406 




Bu. 


mp 


14.0 


SST24 


J143358 


.0+332607 


218 


.4916382 


33 


.4355431 


2.414 


Pt 




r-law 


> 16.5 


SST24 


J143407 


.4+343242 


218 


.5311125 


34 


.5451361 


3.791 




Bu: 




> 15.7 


SST24 


J143410 


.6+332641 


218 


.5445557 


33 


.4447975 


2.263 


Pt 




r-law 


14.1 


SST24 


J143411 


.0+331733 


218 


.5457833 


33 


.2924194 


2.656 


p, 




r-law 


13.8 


SST24 


J143424 


.4+334543 


218 


.6019135 


33 


.7619972 


2.263 


Pt 




r-law 


> 15.2 


SST24 


J143447 


.7+330230 


218 


.6988373 


33 


.0417976 


1.78 


p, 




r-law 


> 17.0 


SST24 


,1143458 


.9+333437 


218 


.7454834 


33 


.5770416 


2. 150 




Bu 




14.2 


SST24 


J143502 


.9+342658 


218 


.7622208 


34 


.4496611 


2.10 




Bu: 


mp 


14.2 


SST24 


J143503 


.2 + 340243 


218 


.7635042 


34 


.0454417 


1.97 




Bu. 


mp 


15.3 


SST24 


,1143504 


.1 + 354743 


218 


.7672272 


35 


.7955055 


2.13 


Pt 




r-law 


16.2 


SST24 


J143508 


.4+334739 


218 


.7854614 


33 


.7942467 


2.10 


P, 




r-law 


15.3 


SST24 


J143509 


.7+340137 


218 


.7904500 


34 


.0269583 


2.080 


Pt 




r-law 


14.6 


SST24 


J143518 


.8+340427 


218 


.8285065 


34 


.0741196 


1.996 




Bu: 


mp 


13.9 


SST24 


,1143520 


.7+340602 


218 


.8361969 


34 


.1007767 


1.730 




Bu. 


mp 


13.8 


SST24 


,1143520 


.7+340418 


218 


.8364868 


34 


.0716324 


1.790 


P> 




r-law 


15.8 


SST24 


J143523 


.9+330706 


218 


.8497772 


33 


.1186829 


2.59 


P. 




r-law 


15.3 


SST24 


J143539 


.3 + 334159 


218 


.9140167 


33 


.6998062 


2.62 


P> 




r-law 


> 16.8 


SST24 


J143545 


.1 + 342831 


218 


.9378204 


34 


.4752998 


2.50 




Bu: 


mp 


> 16.0 


SST24 


J143631 


.8 + 350210 


219 


.1326141 


35 


.0360146 


1.689 




Bu. 


mp 


15.0 


SST24 


J143632 


.7+350515 


219 


.1362610 


35 


.0877495 


1.743 


P. 




r-law 


14.3 


SST24 


J143634 


.3 + 334854 


219 


.1430206 


33 


.8151054 


2.267 


P> 




r-law 


14.9 


SST24 


J143641 


.0+350207 


219 


.1708542 


35 


.0353083 


1.948 




Bu: 


mp 


14.0 


SST24 


J143641 


.6+342752 


219 


.1735382 


34 


.4644394 


2.752 


P> 




r-law 


14.9 


SST24 


J143644 


.2 + 350627 


219 


.1842804 


35 


.1075211 


1.95 


P. 




r-law 


15.6 


SST24 


J143701 


.9+344630 


219 


.2582875 


34 


.7751167 


3.04 




Bu 


mp 


> 15.6 


SST24 


J143725 


.1 + 341502 


219 


.3548889 


34 


.2506104 


2.50 


P. 




r-law 


> 16.2 


SST24 


J143740 


.1 + 341102 


219 


.4176636 


34 


.1841354 


2.197 


P> 




r-law 


14.5 


SST24 


J143742 


.5 + 341424 


219 


.4276276 


34 


.2403145 


1.901 


P. 




r-law 


15.0 


SST24 


J143808 


.3+341016 


219 


.5347443 


34 


.1708908 


2.50 


P> 




r-law 


15.5 


SST24 


J143816 


.6 + 333700 


219 


.5695038 


33 


.6167984 


1.84 




Bu: 


mp 


14.5 



a Redshifts are fro m either Spitzer/IRS (two-decimal point precision; |Houck et"im 120051 ; 
Weed miin et al . 2006a ) or Keck (three decimal point precision ; Soifer ct al . , in prep . J spectroscopy. 
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Optical, Near-IR, and Mid-IR Photometry of DOGs. All fl 



Table 2 

ix densities given in units of fJ.Jy. Note that all measurements and theii 
whether the measurement is statistically significant. 



uncertainties are reported, regardless of 



SST24 


.1142538 


.2 + 351855 


-0 


.05 + 0.04 





.05 + 0. 


10 


SST24 


J142541 


.3+342420 





19 + 


06 





.37+0. 


15 


SST24 


J142554 


.9+341820 





19 + 


.04 





.30 + 0. 


13 


SST24 


J142607 


.8+330425 





.12 + 


.05 





.38 + 0. 


05 


SST24 


J142622 


.0+345249 





44+0 


05 





.52 + 0. 


12 


SST24 


J142626 


.4+344731 





.04+0 


05 





.00 + 0. 


13 


SST24 


J142637 


.3+333025 





.10+0 


05 





.14+0. 


18 


SST24 


J142644 


.3+333051 





.08 + 


.06 





.52 + 


19 


SST24 


J142645 


.7+351901 





.04 + 


.03 





.07 + 


07 


SST24 


J14264S 


.9+332927 





.34+0 


.06 





.52 + 0. 


21 


SST24 


J142652 


.5+345506 





09 + 


.04 





.24+0. 


11 


SST24 


J142653 


.2+330221 





.10 + 


.04 





.18 + 0. 


08 


SST24 


J142724 


.9+350824 





09 + 


04 





.15 + 0. 


13 


SST24 


J 142 748 


.4+344851 


1 


66 + 


06 


1 


26 + 0. 


13 


SST24 


J142759 


.8+351243 





.34+0 


.04 





.40 + 0. 


21 


SST24 


,1142800 


.6+350455 





.40 + 


05 





.51 + 0. 


14 


SST24 


,1142804 


.1 + 332135 





.00 + 


.05 


-0.01 + 0.09 


SST24 


J142810 


.5+352509 





.14+0 


03 


0.32 + 


11 


SST24 


J142814 


.2+352245 





.20 + 


.03 





50 + 0. 


11 


SST24 


J142815 


.4+324720 





33 + 


.04 





.51 + 0. 


07 


SST24 


J142827 


.9+334550 





.20 + 


. 04 


0.22 + 0. 


09 


SST24 


J142832 


.4+340849 





29 + 


02 


0.68 + 0. 


15 


SST24 


,1142842 


.9+342409 


1 


.12 + 


.06 


1 


23 + 0. 


17 


SST24 


J142846 


.6+352701 





.10 + 


.05 





18 + 0. 


12 


SST24 


J142901 


.5+353016 





39 + 


.04 


0.24+0. 


12 


SST24 


J142920 


.1 + 333023 





22 + 


.06 


0.53 + 0. 


09 


SST24 


,1142924 


.8+353320 





08 + 


.04 


0.09 + 0. 


10 


SST24 


J142928 


.5+350841 


-0 


.01 + 0.07 





14+0. 


14 


SST24 


J142931 


.3+321828 


-0 


.13 + 0.07 


-0.09 + 


.12 


SST24 


J142934 


.2+322213 


0.61 + 0.04 


0.39 + 0.06 


SST24 


J142941 


.0+340915 


0.05 + 0.04 


-0 


.07+0 


.16 


SST24 


,1142951 


.1 + 342042 


0.24+0.04 





.26 + 0. 


16 


SST24 


,1142958 


.3+322615 


0.20 + 0.04 





.27+0. 


09 


SST24 


J143001 


.9+334538 


0.28 + 0.06 





.52 + 0. 


12 


SST24 


,1143020 


.4+330344 


0.06 + 0.05 





.17+0. 


09 


SST24 


J143022 


.5+330029 


0.01 + 0.05 





.16 + 0. 


09 


SST24 


J143025 


.7+342957 


0.46 + 0.04 





.70 + 0. 


12 


SST24 


J143028 


.5+343221 


0.35 + 0.06 





.47+0. 


10 


SST24 


J143102 


.2+325152 


-0.04+0.04 





.16 + 0. 


12 


SST24 


J143109 


.7+342802 





02 + 


05 





.24+0. 


09 


SST24 


J143135 


.2+325456 





.41 + 


04 





.80 + 0. 


10 


SST24 


J143137 


.1 + 334501 





.17+0 


06 





.28 + 


12 


SST24 


J143152 


.3+350030 





.14+0 


.03 





.30 + 0. 


09 


SST24 


,1143201 


.8+340408 





43 + 


05 





.43 + 0. 


18 


SST24 


J143216 


.8+335231 





.11 + 


.05 





.23 + 0. 


12 


SST24 


,1143225 


.3+334716 





.07+0 


04 





.13 + 0. 


12 


SST24 


J143242 


.5+342232 





.05 + 


.04 





.20 + 0. 


12 


SST24 


J143251 


.8+333536 





06 + 


05 





.13 + 0. 


13 


SST24 


J143312 


.7+342011 





.53 + 


04 





.57+0. 


13 


SST24 


J143315 


.1 + 335628 





.42 + 


.07 





.73 + 0. 


08 


SST24 


J143318 


.8+332203 





28 + 


05 





.27+0. 


06 


SST24 


,1143321 


.8+342502 





.18 + 


.07 





.50 + 0. 


09 


SST24 


,1143324 


.3+334239 





.24+0 


06 





.44+0. 


12 


SST24 


J143325 


.8+333736 





.20 + 


.06 





.55 + 0. 


10 


SST24 


,1143330 


.0+342234 





.43 + 


05 





.66 + 0. 


12 


SST24 


J143331 


.9+352027 





.18 + 


.03 





.47+0. 


06 


SST24 


,1143332 


.5 + 332230 





09 + 


.05 





.14+0. 


07 


SST24 


J143335 


.9+334716 





36 + 


06 





.39 + 0. 


14 


SST24 


J143349 


.5+334601 





.12 + 


.06 





.26 + 0. 


14 


SST24 


,1143353 
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Table 3 

Median and inter-quartile M* values for PL DOGs, bump DOGs, 
and SMGs using the CB07 library and Chabrier IMF. SFH as noted. 



SFH 


PL DOGs 
Iog(Af,/M s ) 


Bump DOGs 
Iog(Af./Af s ) 


SMGs 
Iog(Af./M s ) 


Instantaneous Burst a 
Major Merger 10 
Major Merger 
Smooth Accretion d 
Smooth Accretion 


in 7i + ' 40 

lu - ' 1 -0.34 

io.9o±g;|g 
ii.06l8:^ 

10.94±g-Jg 

n.2o±°:lg 


10.62±° ;i 
10.741°;^ 
10.88±° ;« 
10.75l°j° 5 

n.03±g:iS 


10-42l° ;f 6 

io.59l8; 3 3 4 6 

10.86+8^ 
10.64ig;iJ 

n-02±8:f 7 



a Simple stellar population 
b SFH from Narayanan ct al. 
c SFH from [Narayanan c tafl ([20lC 
to 2 = 2 (i.e., the peak SFR period) 
d SFH from lDave eFaLl (|20lol) 
c SFH from 
2 ~ 2 — 3 range 



and restricting the time range 



Dave et "all 1 20101) and restricting the time range to the 
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Figure 8. Luminosity per unit frequency as a function of rest-frame wavelength for power-law DO Gs ( red circles). Gray circles indicate 
that Bw and IRAC 5.8^tm and 8. (him data are not used to constrain the SPS models (see section l3.l) . The best-fit CB07 synthesized 
stellar population (assuming a Chabrier IMF and a simple stellar population SFH) is shown in green with the best-fit parameters printed 
in the bottom right of each panel. The inset shows the stellar mass probability density function (marginalizing over model age and Ay) 
for each source. 



DOG & SMG SFHs 



19 





10 34 




0.3 




1Q33 




0.2 


1 






0.1 


M 

X 


10 32 




0.0 








9 


to 


10 31 






CP 








<u 


10 30 






_l 


10 29 




• 




10 28 








0.1 






0.1 



0.1 



Av = 2.4 
Age = 20 Myr 
lbg(M./M Q ) = 1 1.0 
SST24 J I42648. 9+332927 



1.0 

Rest Wavelength (/im) 



1.0 

Rest Wavelength (/im) 



log(M./Mo) 

Av = 1.3 
Age = 160 Myr 
log(M./M ) =11.3 
SST24 J 142800.6+350455 



1.0 

Rest Wavelength (fim) 



A 



log(M./M ) 

Av = 1.8 
Age = 10 Myr 
log(M./Mo) = 10.7 
SST24 J142814.2+352245 




1.0 

Rest Wavelength (/im) 







i — i — i 1 1 — 


0.3 






0.2 






0.1 


ft : 




0.0 




• 


9 


10 11 12 13 • 




log(M./M s ) 








-TT^SAv = 


• 




Age = 10 Myr 




log(M-/M©) =10.1 " 




SST24 J142748.4+344851 
i i ■ i ■ i i i i i i 






10 33 




0.3 








0.2 


i 


10 32 




0.1 










i 


103' 




0.0 

! 


i 

in 








CP 


10 30 






0) 








» 


1029 






_l 










10 28 








0.1 






log(M./M ) 

Av = 2.0 

Age = 30 Myr 
ldg(M./M ) = 10.4 
S$T24 J 142653.2+330221 



1.0 

Rest Wavelength (/im) 



10 33 

n 10 32 

T « 10 3 ' 

Jd 10 30 

J 10 29 

1 2S 1= 
0.1 




Av = 2.0 
Age = 160 Myr 
log(M./M Q ) = 1 1 .3 
SST24 J142759.8+351243 



1.0 

Rest Wavelength (/im) 




■ 

Av = 1.8 
Age - 40 Myr 
log(M./Mo) = 10.7 
SST24 J 1428 10.5+352509 



1.0 

Rest Wavelength (/im) 



10 33 - 
r~ 10 32 - 

I 

N 

* 10 31 - 

I 

(/) 

o. 10 30 - 
t 10 29 - 

_l 

10 28 - 
0.1 




log(M./M ) 

Av = 1.4 
Age - 30 Myr 
log(M./M Q ) = 10.4 
SST24 J142815.4+324720 



1.0 

Rest Wavelength (/im) 



Figure 8. Continued. 
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Figure 8. Continued. 
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Figure 8. Continued. 
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Figure 8. Continued. 
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Figure 8. Continued. 
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Figure 8. Continued. 
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Figure 9. Continued. 
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Figure 9. Continued. 
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Figure 10. Same as Figure |U but for SMGs. 
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Figure 10. Continued. 
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Figure 10. Continued. 
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Figure 10. Continued. 
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Figure 10. Continued. 
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Figure 10. Continued. 
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